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The Casimir force for charge-neutral, perfect conductors of non-planar geometric configurations have been 
investigated. The configurations are: (1) the plate-hemisphere, (2) the hemisphere-hemisphere and (3) the 
spherical shell. The resulting Casimir forces for these physical arrangements have been found to be attractive. 
The repulsive Casimir force found by Boyer for a spherical shell is a special case requiring stringent material 
property of the sphere, as well as the specific boundary conditions for the wave modes inside and outside of the 
sphere. The necessary criteria in detecting Boyer's repulsive Casimir force for a sphere are discussed at the end 
of this investigation. 
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Polder faced a similar problem in studying forces between po- 
larizable neutral molecules in colloidal solutions. Colloidal 
solutions are viscous materials, such as paint, that contain 
micron-sized particles in a liquid matrix. It had been thought 
that forces between such polarizable, neutral molecules were 
governed by the van der Waals interaction. The van der Waals 
interaction is also referred to as the "Lennard-Jones interac- 
tion." It is a long range electrostatic interaction that acts to 
attract two nearby polarizable molecules. Casimir and Polder 
found to their surprise that there existed an attractive force 
which could not be ascribed to the van der Waals theory. Their 
experimental result could not be correctly explained unless 
the retardation effect was included in van der Waals' theory. 
This retarded van der Waals interaction or Lienard-Wiechert 
dipole-dipole interaction is now known as the Casimir-Polder 
interaction QUEUE!!. Casimir, following this first work, 
elaborated on the Casimir-Polder interaction in predicting the 
existence of an attractive force between two electrically neu- 
tral, parallel plates of perfect conductors separated by a small 
gap 1 7]. This alternative derivation of the Casimir force is in 
terms of the difference between the zero-point energy in vac- 
uum and the zero-point energy in the presence of boundaries. 
This force has been confirmed by experiments 1 sf [HJ Hill 
and the phenomenon is what is now known as the "Casimir 
Effect." The force responsible for the attraction of two un- 
charged conducting plates is accordingly termed the "Casimir 
Force." It was shown later that the Casimir force could be 
both attractive or repulsive depending on the geometry and 
the material property of the conductors 11121 U31H4I1 . 



I. INTRODUCTION 
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When two electrically neutral, conducting plates are placed 
parallel to each other, our understanding from classical elec- 
trodynamics tells us that nothing should happen to these 
plates. The plates are assumed to be that made of perfect 
conductors for simplicity. In 1948, H. B. G. Casimir and D. 
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Figure 1 : The vacuum polarization of a photon. 

The Casimir effect is regarded as macroscopic manifesta- 
tion of the retarded van der Waals interaction between un- 
charged polarizable molecules (or atoms). Microscopically, 
the Casimir effect is due to interactions between induced mul- 
tipole moments, where the dipole term is the most dominant 
contributor if it is non- vanishing. Therefore, the dipole inter- 
action is exclusively referred to, unless otherwise explicitly 
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stated, throughout this investigation. The induced dipole mo- 
ments can be qualitatively explained by the concept of "vac- 
uum polarization" in quantum electrodynamics (QED). The 
idea is that a photon, whether real or virtual, has a charged par- 
ticle content. Namely, the internal loop, illustrated in Figure 
3] can be e + e~ , , t + t~ , ir + ir~ or qq pairs, etc. Its cor- 

rectness have been born out from the precision measurements 
of Lamb shift El III and photoproduction O El EEl 
experiments over a vast range of energies. For the almost zero 
energy photons considered in the Casimir effect, these pairs 
last for a time interval At consistent with that given by the 
Heisenberg uncertainty principle AE ■ At = h, where AE is 
the energy imbalance and h is the Planck constant. These vir- 
tual charged particles can induce the requisite polarizability 
on the boundary of the dielectric (or conducting) plates which 
explained the Casimir effect. However, the dipole strength is 
left as a free parameter in the calculations because it cannot 
be readily calculated |21]. Its value can be determined from 
experiments. 



Once this idea is taken for granted, one can then move for- 
ward to calculate the effective, temperature averaged, energy 
due to the dipole-dipole interactions with the time retardation 
effect folded in. The energy between the dielectric (or con- 
ducting) media is obtained from the allowed modes of elec- 
tromagnetic waves determined by the Maxwell equations to- 
gether with the boundary conditions. The Casimir force is 
then obtained by taking the negative gradient of the energy in 
space. This approach, as opposed to full atomistic treatment 
of the dielectrics (or conductors), is justified as long as the 
most significant photon wavelengths determining the interac- 
tion are large when compared with the spacing of the lattice 
points in the media. The effect of all the multiple dipole scat- 
tering by atoms in the dielectric (or conducting) media simply 
enforces the macroscopic reflection laws of electromagnetic 
waves. For instance, in the case of the two parallel plates, 
the most significant wavelengths are those of the order of the 
plate gap distance. When this wavelength is large compared 
with the interatomic distances, the macroscopic electromag- 
netic theory can be used with impunity. The geometric con- 
figuration can introduce significant complications, which is 
the subject matter this study is going to address. 



In order to handle the dipole-dipole interaction Hamilto- 
nian in this case, the classical electromagnetic fields have to 
be quantized into the photon representation first. The photon 
with non-zero occupation number have energies in units of 
hui; where h is the Planck constant divided by 2tt, and u, the 
angular frequency. The lowest energy state of the electromag- 
netic fields has energy Huj/2. They are called the vacuum or 
the zero point energy state, and they play a major role in the 
Casimir effect. Throughout this investigation, the terminol- 
ogy "photon" is used to represent the entity with energy hui, 
or the entity with energy hui/ 2 unless explicitly stated other- 
wise. With this in mind, the quantized field energy is written 
as 



Hn. h — 



hcOi 
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where c is the speed of light in empty space, Q^i is the degree 
of freedom in polarization, k[ is the wave number, n; is the 
wave mode number and Li, the boundary length. The sub- 
script b of Ti' n b denotes the bounded space. For electromag- 
netic waves, the degree of freedom in polarization, 0^/ = 2. 
Similarly, in free space, the field energy is quantized in the 
form 



nil _ n_l/ 



A (ii) A (£2) h (is) 



[k' i {n i ,L i )fdk' 1 dk' 2 dk' 3 , (2) 



where the subscript u of 7i' n u denotes free (unbounded) 
space and the functional fi (Li) in the denominator is equal 
to Czeron^ 1 L^ 1 for a given Li. Here ( zero is the zeroes of the 
function representing the transversal component of the elec- 
tric field. The corresponding zero point energy for bounded 
(or unbounded) space is obtained by setting n s = in equa- 
tions and 0, respectively. 



II. REFLECTION DYNAMICS 

In principle, the atomistic approach utilizing the Casimir- 
Polder interaction |2] explains the Casimir effect observed in 
any system. Unfortunately, the pairwise summation of the in- 
termolecular forces for systems containing large number of 
atoms can become very complicated. H. B. G. Casimir, real- 
izing the linear relationship between the field and the polariza- 
tion, devised an easier approach to the calculation of Casimir 
effect for large systems such as two perfectly conducting par- 
allel plates 1 7]. The Casimir effect have been also explained 
by Schwinger utilizing his original invention, "Source The- 
ory" 

In this investigation, we do not follow Casimir's energy 
method, nor do we take the route of Schwinger's source the- 
ory. Instead, we adopt the vacuum pressure approach intro- 
duced by Milonni, Cook and Goggin 1 22], which is a simple 
elaboration on Casimir's original calculation technique utiliz- 
ing the boundary conditions. We choose to consider the vac- 
uum pressure approach over both Casimir's energy method 
and Schwinger's source theory not because it is a superior 
technique, but simply because it is the easiest one for the phys- 
ical arrangements considered in this investigation. 

The three physical arrangements for the boundary config- 
urations considered in this investigation are: (1) the plate- 
hemisphere, (2) the hemisphere-hemisphere and (3) a sphere 
formed by brining two hemispheres together. Because the ge- 
ometric configurations of items (2) and (3) are special versions 
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Figure 2: The plane of incidence view of plate-hemisphere configu- 
ration. The waves that are allowed through internal reflections in the 

hemisphere cavity must satisfy the relation A < 2 -R'2 — R'i ■ 



of the more general, plate-hemisphere configuration, the ba- 
sic reflection dynamics needed for the plate-hemisphere case 
is worked out first. The results can then be applied to the 
hemisphere-hemisphere and the sphere configurations later. 

The vacuum-fields are subject to the appropriate bound- 
ary conditions. For boundaries made of perfect conduc- 
tors, the transverse components of the electric field are zero 
at the surface. For this simplification, the skin depth of 
penetration is considered to be zero. The plate-hemisphere 
under consideration is shown in Figure The solutions 
to the vacuum-fields are that of the Cartesian version of 
the free Maxwell field vector potential differential equation 

V 2 A (Rj - c~ 2 d}A (r) = 0, where the Coulomb gauge 
V • A = and the absence of the source p for scalar potential 



$ have been imposed, $ ( p 



R 



0. The electric and the 



magnetic field component of the vacuum-field are given by 
E = —cr x dtA and B = V X A, where A is the free field vec- 
tor potential. The vanishing of the transversal component of 
the electric field at the perfect conductor surface implies that 

the solution for E is in the form of E oc sin ( 2ttX~ 



where A is the wavelength and 



L 



is the path length between 



the boundaries. The wavelength is restricted by the condition 



A < 2 



R' 



R'l 



2^2, where R'2 and R\ are two im- 
mediate reflection points in the hemisphere cavity of Figure|2] 
In order to compute the modes allowed inside the hemisphere 
resonator, a detailed knowledge of the reflections occurring 
in the hemisphere cavity is needed. This is described in the 
following section. 



A. Reflection Points on the Surface of a Resonator 

The wave vector directed along an arbitrary direction in 
Cartesian coordinates is written as 



(3) 



where 



fc'i (fci,x! k'i,y> fci,z) — fcl.i^i, 

i=l 

{i = l—> k' lx , e"i = x, 
i = 2^ k' ly , e 2 = y, 
i = 3 — > k' l z , e 3 = z. 

- -1 3 

Hence, the unit wave vector, k 1 1 — fc'i 2i=i M i^i- De- 
fine the initial position R'o for the incident wave fc'i, 
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(4) 



i=i 



where 
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i = 3 



0,rc' 
r 0,y> 
r 0,z- 



Here it should be noted that R'q really has only components 
r' Q x and r' Q z . But nevertheless, one can always set r' = 
whenever needed. Since no particular wave vectors with 
specified wave lengths are prescribed initially, it is desirable 
to employ a parameterization scheme to represent these wave 
vectors. The line segment traced out by this wave vector k\ 
is formulated in the parametric form 



r'i = eifc'i + R'o 

3 







-1 




fc'i 


k' 

K t,i 



(5) 



where the variable ^1 is a positive definite parameter. Here 
R'l is the first reflection point on the hemisphere. In terms of 
spherical coordinate variables, R\ takes the form 



R'i {r'i, 0'i, 4>i) = r 'iYl A 'i,i e "i' 



(6) 



i=l 



where 



l) 



A j 1 = sin 0[ cos 
A^ 2 = sin 9[ sin (/>[ , 
A' 13 = cos9[. 
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Here r\ is the hemisphere radius, 6[ and (f>[ are the polar and 
the azimuthal angle respectively of R'i at the first reflection 
point. Notice that subscript i of r[ denotes "inner radius" not 
a summation index. 

By combining equations (|5} and (|6j, we can solve for the 
parameter £1. It can be shown that 



6 = £i* 



(7) 



where the positive root for £i have been chosen due to the 
restriction £i > 0. The detailed proof of equation |7} is given 
in Appendix A of reference Qj. Substituting £i in equation 
(|5}, the first reflection point off the inner hemisphere surface 
is expressed as 



R 



r 'o,i + 



k'i 



h' 

K l,i 



(8) 



where £i. p is from equation Q. 

The incoming wave vector k\ can always be decomposed 
into parallel and perpendicular components, fe'^ii and k'i t ±, 
with respect to the local reflection surface. It is shown in Ap- 
pendix A of reference |1] that the reflected wave vector k' r 
has the form 



k' 



n' x k'i 



x n' — a r it n' ■ k'in', 



where the quantities a r n and a r> ± are the reflection coeffi- 
cients and n' is a unit surface normal. For the perfect reflect- 
ing surfaces, a r _\\ = a r ,± — 1. In component form, 

3 

k' r = {"n-L [n' n k' i in ' n - n[k' hn n' n ] 
i=i 

- a r,\\ n 'n k i,n n l} e "'> 

where it is understood that n' is already normalized and Ein- 
stein summation convention is applied to the index n. The sec- 
ond reflection point R'2 is found then by repeating the steps 
done for R\ and by using the expression k' r = k' r / k'. 

R'2 = R'l + &,pk'r 

- R'i + 6, P 



a r ,± 


n' x k'i 


x n' 


— o. r \\n' 


■ k'in' 


a r ,± 


n' x k\ 


x n' 


— oi r \\n' 


■ k'in' 



where £ 2 ,p is me new positive definite parameter for the sec- 
ond reflection point. 

The incidence plane of reflection is determined solely by 
the incident wave k'i and the local normal n'i of the reflect- 
ing surface. It is important to recognize the fact that the sub- 
sequent successive reflections of this incoming wave will be 
confined to this particular incidence plane. This incident plane 



can be characterized by a unit normal vector. For the system 
shown in Figure|5J 



and 



k % — k'i, 
n'n',i = -£i,pfe'l - -R'o- 



The unit vector which represents the incidence plane is given 
by 



n'pA 



n ' p ' 1 | t ijkki jr o k €j 1 

i=l 



where the summations over indices j and k are implicit. If 
the plane of incidence is represented by a scalar function 
/ (x' , y' , z') , then its unit normal vector n' Pt i will satisfy the 

relationship n' p i oc V'/p.i ( x \ y\ z ') ■ It can be shown (see 
Appendix A of reference 1 1]) that 



fp.l Ol,^!^) 

where 



(9) 



2=1 




and 



-00 < {v[ = x', v' 2 = y', v' z = z'} < 



The surface of a sphere or hemisphere is defined through 
the relation 

f hemi (x',y',z')-K] 2 -EM 2 ' 

i=l 

where r\ is the radius of sphere and the subscript i denotes 
the inner surface. The intercept of interest is shown in Figure 
[3] The intersection between the hemisphere surface and the 
incidence plane f p i (v[, v' 2l ^3) is given through the relation 

hemi ( x \ y', z ') - fp,l ( x ', y\ z ') = 0. 

After substitution of f p . x (x', y', z') and f heml (x' , y' , z') , we 
have 



i=l 

where 







{k'] 2 - 


n' p> i 



n' p ,i eijkki jr' k4 } - K 
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Intercept 




The x coordinate is coming out of the page! 



Figure 3: The thick line shown here represents the intersection be- 
tween hemisphere surface and the plane of incidence. The unit 
vector normal to the plane of incidence is given by n' p ,i = 



the initial reflection point R' x can be expressed in terms of the 
spherical coordinate variables {r' i: 9' x , </>[) , 

3 
i=l 

where 

1 — ► v' x x = r\ sin 9[ cos <f>' x , 

3 -> *4, 3 = r i cos6l 'i- 

Here r[ is the hemisphere radius, 4>' x and (9' l7 the polar and 
azimuthal angle, respectively. The angles (j>[ and 9[ are found 
to be (see Appendix A of reference 1 1 ]): 

lim e ^o (0 < 4>' x <±ir- |e|) , 
lim £ ^o (tt + |e| < < §tt - |e|) , 



6^ = arctan 



'r?,-T ll3 ||n'p 



i||™p.i|| £iikk ld r ok 



(12) 



The term [rj] can be rewritten in the form 



where 



M] 9 = £[<<] a . 





= 1 


'<! 


= r' 

2, a; 


< <2 


= r' 

1,3/ 


> r i,3 


= r' 



Solving for i/, it can be shown (see Appendix A of reference 
CI) that 



n p.i 



± 



"1 




2 





(10) 



where i = 1,2,3 and e^fc is the Levi-Civita symbol. The 
result for i/ shown above provide a set of discrete reflection 
points found by the intercept between the hemisphere and the 
plane of incidence. Using spherical coordinate representations 
for the variables r\ x , r[ 2 and r\ 3 , 



r 'i l = r i sin (9' cose//, 



r i.2 = »** sin 6' sin 



L r i,3 = r i C0S( 



lim e ^o (|tt + |e| < (j)' x < ir - e) , 
lim £ ^ (|7r+|e| < <j>' x < 2ir - |e|) , 



arctan 



T i,l- T l.l|K'j>.l|| e ijkk' 1:j r' Qtk 



9' x = arctan 



lim^o (0 < 6' x < ivr - |e|) , 



Ti,ieijfcfci.j^o,fc+Ti,2e2mnfci, m r 0| „ 
||n'j>,i||{ T i,3- T i.3|K'j>,i|r 1 e3 Se rfci, 9 rJ jr ] 

lim e ^o + \s\ < 9[ <n), 
9[ = arctan 



t Ll±IL2 



where 



k x i — 1,2,3. 



(13) 



(14) 



T 1,3- T 1,3||«'p,i|| ^rM.^o,, (15) 
T l,l £ ljfc fc l,j r 0,fc+ T l,2e2m ra fci, m r 0| „ 

||™'p.i||{ T f,3~ T i,3||n'p,i|| 1 e 3 qrk' 1:q r' r j 



Similarly, the second reflection point on the inner hemisphere 
surface is given by (see Appendix A of reference 1 1J): 

3 

R'2 (r'i,9' 2 ,(j)' 2 ) = ^ V 2,i ( r 'i> 02, 02 ) e ~i; ( i6 ) 



where 



l-n/ 21 = r<sm0 2 cos0 2 , 



'2,3 



r; sin t/, sm ( 



r'i cos V' 2 . 



Here the spherical angles cf>' 2 and 9' 2 are found to be (see Ap- 
pendix A of reference 01): 
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# 2 = arctan 



n'p.i 



n' P A 



€ ljkk' 1 ,j r 0,k 



n' Pl i e 2mnk' lm r' 0n i> 22 



1 




2 


n 'p,i 



1/2 



e 3qrk[ q r' Q r j> 2 ,3 



e 3qrk' lq r' , 



-1/2N 



lilHe^O (0 < 2 < §7r - | E |) , 

lim £ ^o (7T+|e| < 02 < |tt — |er|) , 



/>' 2 = arctan 



1 




2 





e 2mnk' lm r' 0n i> 22 



n' Pi i 



e 2mn^ ,m r 0,n 



e iife^ij r d,fc ^2,1 



-l/2\ 



(17) 



lim e _ (i7r+|e| < 0' < tt) 



9 2 = arctan — 





1 






2 




1 




-l 



i 



eijkkij' 1 



0,k 



1 




2 


«'p,l 



1 




2 





e 2mnk[ m r' n i> 2 2 



e 3grfci i „''o,r ^2,3 



-1/2N 



71 P,l 



e 3qrk' lq r' Qr 



where the variables have the definition: 
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where k[ x — k' x , , k[ 2 = k' y , and fc^ 3 = k' z , ; 



and 



r li2 = K] 2 sin (tt - 26 mc ) , 
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It can be shown that Nth reflection point inside the hemi- 
sphere (see Appendix A of reference | j[) is, 

3 

R' N (rj, 6' N , <j>' N ) = V v' N>i {r'^e'x, <j>' N ) e h 



where 
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(25) 



The angular variables 6*^, and 0jy corresponding to iVth re- 
flection point ii'jv are given (see Appendix A of reference 
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Figure 4: The surface of the hemisphere-hemisphere configuration 
can be described relative to the system origin through R, or relative 
to the hemisphere centers through R'. 

cases, the preferred choice for the system origin, from which 
the variables are defined, depend on the physical arrangements 
of the system being considered. For a sphere, the natural 
choice for the origin is its center from which the spherical 
variables ,cf>') are prescribed. For more complicated 

configuration shown in Figure |4j the preferred choice for ori- 
gin really depends on the problem at hand. For this reason, 
a set of transformation rules between (r^, 8', 0') and (r^, 6, (f>) 
is sought. Here the primed set is defined relative to the sphere 
center and the unprimed set is defined relative to the origin of 
the global configuration. In terms of the Cartesian variables, 
the two vectors R and R' describing an identical point on the 
hemisphere surface are expressed by 



(30) 



-l/2\ 



where 



(29) 

where {v'ma ■ i = 1,2, 3} are found by simply replacing in 



equations \22\ . d23l and J24> the following: 

Ci (IY 2 ) - Ci (Ti,jv) , 

C2 (r 1>2 ) -> c 2 (T hN ) , 

I Ca (Ti, 2 ) -> C 3 (Ti,n) , 
where IYjv is 

r 1 , 2V = H] 2 sin([AT-l] [n~29 mc ]). 

The details of all the work shown up to this point can be found 
in Appendix A of reference (J. 

The previously shown reflection points (R'i, R'2 and R'n) 
were described relative to the hemisphere center. In many 



R(f U 1/2,1^3) = E i= i^ei, 

(1^1, f 2, ^3) (x,y, z) , 

(ei,e 2 ,e3) -> (x,y, z) . 

The vectors R and Rl are connected through the relation 

3 

R{vy,V2,vz) [i^r,j + fi] e A j 

with 



representing the position of hemisphere center relative to the 
system origin. One has then the condition 



\v% - v T ^ - i/f] Si = 0. 
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In terms of spherical coordinate representation for (u\, 1/2, 1/3) 
and (u[, u' 2 , i/' 3 ) , we can solve for 9 and <\> to yield the results 
(see Appendix B of reference iQ]): 



arctan 



Vt 2 + t'a sin 6' sin 0' 



r- sin y cos < 



(31) 



r i5 ,<j) ,Rt 
arctan 



i*r 1 + v t.2 + sin 9' [cos 0' + sin < 



2*T,3 + H COS 9' 



x < cos arctan 



sin arctan 



/ ^t,2 + r i sm 0' sin <// 
\pt 1 + H sin 6*' cos 0' 



CT2 + H sin 0' sin 0' 



, (32) 



where the notation and indicates that and 9 are explicitly 
expressed in terms of the primed variables, respectively. It is 
to be noticed that for the configuration shown in Figure 0] 
the hemisphere center is only shifted along y by an amount 
of vt,2 = a, which leads to = 0. Nevertheless, the 

derivation have been done for the case where vta ^ 0, % = 
1,2,3 for the purpose of generalization. With the magnitude 

R defined as 



R 



\ i=i 



where 



A; (6', ft) = sin 6' cos </>', 
A' 2 {9',<t>') = sin 0' sin 9/ , 
A' 3 {6') = cos6', 

the vector R ( r', A, A', Rt J is written as 



R 



6 J 

(r> i7 K, A', Rt) = x J2 ^ + r 'i A tf E A * e - < 33 > 

\ i=l i=l 



Here A,, i = 1, 2, 3 is defined as 



= sin f cos ( 



A 2 



= sm f sm ( 



A 3 = cos ( 



The details of this section can be found in Appendices A and 
B of reference {]]] . 



B. Selected Configurations 

Having found all of the wave reflection points in the hemi- 
sphere resonator, the net momentum imparted on both the 
inner and outer surfaces by the incident wave is computed 
for three configurations: (1) the sphere, (2) the hemisphere- 
hemisphere and (3) the plate-hemisphere. The surface ele- 
ment that is being impinged upon by an incident wave would 
experience the net momentum change in an amount pro- 
portional to Ak' i nn er Cj R's,i, R's.o^j on the inner side, and 

Ak' outer ( ) R's.i + a>R's,i \ on the outer side of the surface. 

The quantities Ak' inner an d Ak' ou t er are due to the contribu- 
tion from a single mode of wave traveling in particular direc- 
tion. The notation ( ; R' s ,i, R'sfi^j of Ak' inner denotes that 
it is defined in terms of the initial reflection point i£' s ,i on 

the surface and the initial crossing point R' s ,o of the hemi- 
sphere opening (or the sphere cross-section). The notation 

( ; R's.i + a>R's,i^ °f Ak' ou ter implies the outer surface re- 
flection point. The total resultant imparted momentum on the 
hemisphere or sphere is found by summing over all modes of 
wave, over all directions. 



1. Hollow Spherical Shell 

A sphere formed by bringing in two hemispheres together 
is shown in Figure [5] The resultant change in wave vector 
direction upon reflection at the inner surface of the sphere is 
found to be (see Appendix CI of reference |l|]), 



Ak\ nner R' s i, R' S Q*j 



Ann cos 0,-. 



R 



, 2 (r^ >2 ) -^,1(^,1^) 



R' s ,i, (34) 



where 




Here 9i nc is from equation Mil and, the two vectors 
R s ^i {r l i ,A' s ^j and-R S! 2 {^[1^2) have the generic form 



■5 



(35) 



where 



' K,N,1 (K,N> ( 


K,n) 


= smt 


's,N 




K,n) 


= sin 


9' 

's,N 


K,N,3 ( 


0's,n) 


= cos 


9' 

7 s,N 



COS ( 



sm 



s,N' 
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Figure 5: Inside the cavity, an incident wave k'i on first impact point 
R'i induces a series of reflections that propagate throughout the en- 
tire inner cavity. Similarly, a wave k'i incident on the impact point 
R'i + aR'i, where a is the thickness of the sphere, induces reflected 
wave of magnitude 1 1 fc'i 1 1 . The resultant wave direction in the exter- 
nal region is along R'i and the resultant wave direction in the res- 
onator is along —R'i due to the fact there is exactly another wave 
vector traveling in opposite direction in both regions. In both cases, 
the reflected and incident waves have equal magnitude due to the fact 
that the sphere is assumed to be a perfect conductor. 



The label s has been attached to denote the sphere, and the ob- 
vious index changes in the spherical variables 9' s N and <fr' s N 
are understood from the set of equations (1261 . d27l . J28t and 
(I29> . Similarly, the resultant change in wave vector direction 
upon reflection at the outer surface of the sphere is shown to 
be (see Appendix CI of reference jl||), 



Ak' ou ter (;R's,i+aR' sA ) = 4 



cos 6 inc R' s,!, 

(36) 



where 




The details of this section can be found in Appendix CI of 
reference 1 1 1. 



2. Hemisphere-Hemisphere 

For the hemisphere, the changes in wave vector directions 
after the reflection at a point R'h,i inside the resonator, or af- 
ter the reflection at location R'h.i + aR'h,i outside the hemi- 
sphere, can be found from equations (I34> and (1361 with obvi- 



ous subscript changes, 



Ah', 



H R'h,l, R'h.O^j 

Anir cos 9i nc 



and 



where 



Rh,2 faKQ-Rh,!^',^) 
- (; R'h 



R' 



(37) 



aR' h . 



cos 6 inc R' h A, 



(38) 



< e mc < tt/2, 
n = 1,2, ■ • ■ . 

The reflection location Rh.N {r'i, A-h,Ni A^ N , Rr,h ) has the 
generic form (see Appendix C2 of reference Ql), 

Rh ; jv (r'i, A^jv, A'^jy, Rtm ) 

2 3 



[ VT < h > 



(39) 



where 



Afc,JV,l (®h,N, <frh,N^J 
Aft,AT,2 (&h,N, 4>h,N^ 

A 



sin 6^ cos <ph,N, 
■ sin 9 h<N sin <j) h N , 



kh,N,3 (&h,N^ — COS0h,N- 



The subscript h here denotes the hemisphere. The expressions 
for A' h N j, i = 1,2,3, are defined identically in form. The an- 
gular variables in spherical coordinates, 9h,N an d 4>h,N, can 
be obtained from equations d3 1 1 and i32\ , where the obvious 
notational changes are understood. The implicit angular vari- 
ables, 9' h N and <j>' h N , are the sets defined in equations (1261 
and (I27i for 9' s N , and the sets from equations d28i and ( I29i 
for<p' sN . 

Unlike the sphere situation, the initial wave could eventu- 
ally escape the hemisphere resonator after some maximum 
number of reflections. It is shown (see Appendix C2 of 
referencel 1]) that this maximum number for internal reflec- 
tion is given by 



(40) 



where the notation fZ;, 



'-*h,max\ G denotes the greatest integer 
function, and max is given by 



1 



29i. 



it — arccos 









R'o 















R'o 






R'o 




-'&,})] 



(41) 
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Intercept 




Figure 6: The dashed line vectors represent the situation where only 
single internal reflection occurs. The dark line vectors represent the 
situation where multiple internal reflections occur. 



Here £i. p is from equation Q and 0j nc is from equation i2l\ . 
The above results of Ak' inner ^; R'h,i, R'h,o) an d 

Ak' ou ter R'h.i + aR'h.ij have been derived based on the 
fact that there are multiple internal reflections. For a sphere, 
the multiple internal reflections are inherent. However, for a 
hemisphere, it is not necessarily true that all incoming waves 
would result in multiple internal reflections. Naturally, the 
criteria for multiple internal reflections are in order. If the ini- 
tial direction of the incoming wave vector, k'i, is given, the 
internal reflections can be either single or multiple depend- 
ing upon the location of the entry point in the cavity, R'q. As 
shown in Figure |6j these are two reflection dynamics where 
the dashed vectors represent the single reflection case and the 
non-dashed vectors represent multiple reflections case. Be- 
cause the whole process occurs in the same plane of incidence, 
the vector R' / = —AqR'q where Ao > 0. The multiple or sin- 
gle internal reflection criteria can be summarized by the rela- 
tion (see Appendix C2 of reference 1 1]): 



3 3 r „ 

IN 

[j=i i=i l 

XX. 

U=l 





i 




R'f 


~~ 2 


R'o\\ 




(42) 



Because the hemisphere opening has a radius r\ , the following 




Figure 7: The orientation of a disk is given through the surface unit 
normal n' p . The disk is spanned by the two unit vectors 6' v and cf>' . 



criteria are concluded: 



R' 



<r 



single internal reflection, 



(43) 



> r[, multiple internal reflections, 



where 



is defined in equation J42i . The details of this 
section can be found in Appendix C2 of reference 1 1]. 

3. Plate-Hemisphere 

A surface is represented by a unit vector n' p , which is nor- 
mal to the surface locally. For the circular plate shown in 

Figure^] its orthonormal triad (n' p , §' p , <j)' p ^j has the form 

3 3 BA' 

z— 1 z— 1 P 



1 dA' 



where 



A p,i = sin °'p cos 4>' p , 



V,2 \ u p>Vp) 



sin6>pSiri(/)p, 



(6' p ) = cos9' p . 

For the plate-hemisphere configuration shown in Figure [8] 
it can be shown that the element R p on the plane and its veloc- 
ity dR p /dt are given by (see Appendix C3 of reference 1 1]): 



Rp [A-p> Rr,pj 



, , 9A' p , 

VT,p,i + V p fit 



<<t>' p dA' p , 



-I 2 



80. 



p 



1/2 



XX-< 



(44) 
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V VK ftftl ain ft' f)A)' 






Point on plate 




1/2 


Plate Origin - — i 


















A 





_ d6 p dcb' p p d^ p 



where 



bp = p (6»p, 0p, vt,p,i, vt, p ,2) 

( V T,p,2 



— arctan 



\vt. 



p,i 



Sill COS 0p 



(46) 



# p = bp (o' p , <p' p , Rt,p^ 

v t, p ,i + vt, p ,2 + sin8' p [cos0p + sin</> p ] 



arctan 



cos arctan 



vt, p ,2 + sin 0' p 



cos ( 



sm < 



* V T,p,l 
( VT,p,2 + 



sm 0' cos (p' p 



+ sin arctan , 

\ V V T,P,1 + sm Vp cos Vp 



(47) 



and 



A p ,i (j)p,(/)p^ = sinBp cos 



< A 



p,2 



(48) 



A p ,3 f#p) 



cos 6>„ 



The subscript p of 4> p and # p indicates that these are spherical 
variables for the points on the plate of Figure|8j not that of the 
hemisphere. It is also understood that A' p 3 and A p ,3 are inde- 
pendent of cj) p and 4> p , respectively. Therefore, their differen- 
tiation with respect to 4>' p and (b p respectively vanishes. The 



-V" 



Figure 8: The plate-hemisphere configuration. 



quantities 0' p and cj>' are the angular frequencies, and vt,p,i 
is the translation speed of the plate relative to the system ori- 
gin. The quantities v' P ,e' p and v 1 ' Pj< y are the lattice vibrations 

along the directions Q' p and <f>' p respectively. For the static 



plate without lattice vibrations, v' pA 



and v'p.^i vanishes. 



In the cross-sectional view of the plate-hemisphere system 
shown in Figure [9] the initial wave vector k\ traveling to- 
ward the hemisphere would go through a series of reflections 
according to the law of reflection and finally exit the cavity. 
It would then continue toward the plate, and depending on 
the orientation of plate at the time of impact, the wave-vector, 
now reflecting off the plate, would either escape to infinity or 
re-enter the hemisphere. The process repeats successively. In 
order to determine whether the wave that just escaped from 
the hemisphere cavity can reflect back from the plate and re- 
enter the hemisphere or escape to infinity, the exact location 
of reflection on the plate must be known. This reflection point 
on the plate is found to be (see Appendix C3 of reference 1 1 ]), 



2^i=i a<b' n p,' 



IK'p,i|| 



d K.l ( A ' I £ '"" fc l,m T 0,r. 
2^1 = 1 d6> A I>J ' 



dk' dk' 

del dd>L 



2} 1/2 



h'p.i 



A^A ri B r B 



CpC 1 



(49) 



where the translation parameter vx.p.j — 0; and the terms 

(A c ,B c ,C t: ) , (A 7 ,S 7 ,C 7 ), (Ap,Bp,Cp) and j D are de- 
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fined as (see Appendix C3 of reference 1 1]): 



v-3 9A;,i 

x 1 2^=i ae' 



Ay - £ijkRh,N hima!c ,3 [Rh,N h , maa ,+l,k ~ Rh,N h , 

1 



B~, 



Cy - e ijk [Rh,N himax ,j - Rh,N h>max +l,j] 

x [Rh.N himax +i,k - Rh,N h , max ,k] ; 



Ap = Eli 



AfA 



B I3 - Ei=l ^h,N htrl 



V 3 



1/2 



1 'pM\ E ei 3 kk 'l,j r 0,k C /3 1 C 1 1 



7^/3 „ 



3 r^A' 



9A' 



y^e»jfcAp,jfcjv hi „ iax +i„ 



(=i 

3 



/3°7 



<9A' dA' 

-inr- c c 



m 



V 



i=l 



(50) 



(51) 



(52) 



(53) 



It is to be noticed that for a situation where the translation 
parameter vr,p,j — 0, the A becomes identical to A' in form. 
Results for A can be obtained from A' by a simple replacement 
of primed variables with the unprimed ones. It can be shown 
that the criterion whether the wave reflecting off the plate at 
location R p can re-enter the hemisphere cavity or simply es- 
cape to infinity is found from the relation (see Appendix C3 



Plane of Incidence 




Figure 9: The intersection between oscillating plate, hemi- 
sphere and the plane of incidence whose normal is n' p ,i = 



of reference Ql), 



E 



dA' OA' 

p,s ^ ls -p,s 



2^z=l dd>> \ 2 ^P,i 



d0' p 
2^ 1/2 



3 d K,i 



X 



||n' p ,i I 

AyA g jpRyAp B C B B 



X E i^ 1 - [ k N h , m a, + l^ n 'p,k n p,k 
\k=l 

- n 'p,z k N himax +i.kn' pJ: ] -a r> \\n' Ptk 
xk N h , max +i,kn' p A) , 



C, 
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A 



(54) 



where i = 1, 2, 3 and is the component of the scale vector 
= ^2 i=1 £,n,i£i explained in the Appendix C3 of reference 

m. 

In the above re-entry criteria, it should be noticed that 
Rq < r i ■ This implies r' i < r^, where r\ is the radius of 
hemisphere. It is then concluded that all waves re-entering the 
hemisphere cavity would satisfy the condition = £ Kj 2 = 
£ K .3. On the other hand, those waves that escapes to infin- 
ity cannot have all three £ K ^ equal to a single constant. The 
re-entry condition £ K) i = £ K ,2 = £k,3 is just another way 
of stating the existence of a parametric line along the vector 
k r ,N h max +i that happens to pierce through the hemisphere 
opening. In case such a line does not exist, the initial wave 
direction has to be rotated into a new direction such that there 
is a parametric line that pierces through the hemisphere open- 
ing. That is why all three £ K ; values cannot be equal to a 
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Our formulation of dynamical Casimir effect here, however, 
has no resemblance to that of Schwinger's work to the best of 
our knowledge. This formulation of dynamical Casimir force 
is briefly presented in the following sections. The details of 
derivations pertaining to the dynamical Casimir force can be 
found in Appendix D of reference Qt • 



Figure 10: Because there are more vacuum-field modes in the exter- 
nal regions, the two charge-neutral conducting plates are accelerated 
inward till the two finally stick. 



single constant. The re-entry criteria are summarized here for 
bookkeeping purpose: 

£k,i = £,k,2 = £«,3 '■ wave reenters hemisphere, 

else : wave escapes to infinity, 

(55) 

where else is the case where = £ Kj 2 = 6re,3 cannot be 
satisfied. The details of this section can be found in Appendix 
C3 of reference (l|]. 



1. Formalism of Zero-Point Energy and its Force 

For massless fields, the energy-momentum relation is 

H'„ s = E Tota i = pc, (56) 

where p is the momentum, c is the speed of light, and H' n 
is the quantized field energy for the harmonic fields of equa- 
tion Q for the bounded space, or equation 10 for the free 
space. For the bounded space, the quantized field energy 
H' nB = H' n b of equation Q is a function of the wave num- 
ber k[ (ni) , which in turn is a function of the wave mode value 
n.- L and the boundary functional fi (Li) , where Li is the gap 
distance in the direction of 

Li = R' 2 ■ h - R'i ■ Bi e,. 



C. Dynamical Casimir Force 

The phenomenon of Casimir effect is inherently a dynami- 
cal effect due to the fact that it involves radiation, rather than 
static fields. One of our original objectives in studying the 
Casimir effect was to investigate the physical implications of 
vacuum-fields on movable boundaries. Consider the two par- 
allel plates configuration of charge-neutral, perfect conductors 
shown in Figure 05| Because there are more wave modes in 
the outer region of the parallel plate resonator, two loosely re- 
strained (or unfixed in position) plates will accelerate inward 
until they finally meet. The energy conservation would require 
that the energy initially confined in the resonator when the two 
plates were separated be transformed into the heat energy that 
acts to raise the temperatures of the two plates. 

Davies in 1975 followed by Unruh in 1976 0, have 
asked the similar question and came to a conclusion that when 
an observer is moving with a constant acceleration in vacuum, 
the observer perceives himself to be immersed in a thermal 
bath at the temperature 



T = 



HR 
2-Kck 1 ' 



where R is the acceleration of the observer and k' , the wave 
number. The details of the Unruh-Davies effect can also be 
found in the reference @]. The other work that dealt with 
the concept of dynamical Casimir effect is due to Schwinger 
in his proposals 0, l25ll to explain the phenomenon of sono- 
luminescense. Sonoluminescense is a phenomenon in which 
when a small air bubble filled with noble gas is under a strong 
acoustic-field pressure, the bubble will emit an intense flash 
of light in the optical range. 



Here R\ and R'2 are the position vectors for the involved 
boundaries. As an illustration with the two plate configuration 
shown in Figure0)| R\ may represent the plate positioned at 
z = and R'2 may correspond to the plate at the position 
z = d. When the position of these boundaries are changing in 
time, the quantized field energy TL' n = TL' n b will be modi- 
fied accordingly because the wave number functional k[ (n.i) 
is varying in time, 

dfcj _ 9fc[dnj dfj dL t 

~dt ~ dn~^t Sl [ l) + ni dL,^T 



* it ^ dk 'i ■ 
fi (Li) —^ni 

oni 



dfi f 
ni dLi l 



Here the term proportional to hi refers to the case where the 
boundaries remain fixed throughout all times but the number 
of wave modes in the resonator are being driven by some ac- 
tive external influence. The term proportional to Li represents 
the changes in the number of wave modes due to the moving 
boundaries. 

For an isolated system, there are no external influences, 
hence hi = 0. The dynamical force arising from the time 
variation of the boundaries is found to be (see Appendix Dl 
of reference 1 1||), 




d 2 H' n 

C <*,5 "2 + 0- _ S i°<) 

d[k'i\ 





1" 






1" 




n s - 






n s - 








2. 






2 





■/] (1 - Sij)C at5 nj 



3=1 



dfi d 2 K s ■ 



(57) 
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where C a .i, C a . 2 , C a ^, C aA , C a<5 , C afi and C aJ are de- 
fined as 



r — x^ 3 ^2<l 

_ 2jj=l dk' 



G a .2 = Eli (1 - ([«« + 5] fe ) 2 fc ^ (58) 

I ^, 3 = ELi(i-^)(K + 5] fi ) 2 M 2 ; 



V 



-1/2 



C Q ,5 — 



Cq,iC q ,4 [C^ ^Cq^ Cjj] 



2 [n s + |] ^ 2 C Qi iC„ j2 Cq,4 



^-([n. + l]^) 2 
2 [n, + |] W a . x C a ^ 

C Qj iC Q .3C Q .4 



+ 



C a ,7 



Ca,2Ca,4 

r (K+i] fe ) 2 
c2.i -([». + £] ^) 2 ' 

= CI, - ([n s + ±] hcf 



(59) 



(60) 



(61) 



(62) 



The force shown in the above expression vanishes for the one 
dimensional case. This is an expected result. To understand 
why the force vanishes, we have to refer to the starting point 
equation, 



3 

E 

1=1 







1" 




[( 


n s - 










h 2_ 





-1 awi 



-1/2 



EM 2 



dk> 
dp', = 0, 



(63) 
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Figure 1 1: A one dimensional driven parallel plates configuration. 



of the Appendix Dl of reference (jl- The summation 
here obviously runs only once to arrive at the expression, 
dTL' ris /dk' i = \n s + I] fie- This is a classic situation where 
the problem has been over specified. For the 3D case, equation 
(I63> is a combination of two constraints, Y^=i Wi? an< ^ • 
For the one dimensional case, there is only one constraint, 
H' n . Therefore, equation ( I63> becomes an over specification. 
In order to avoid the problem caused by over specifications 
in this formulation, the one dimensional force expression can 
be obtained directly by differentiating the equation (I56> in- 
stead of using the above formulation for the three dimensional 
case. The ID dynamical force expression for an isolated, non- 
driven systems then becomes 



zt n df OH' ■ . 
c dL dk' 



(64) 



where T' is an one dimensional force. Here the subscript a 
of J 7 ' a have been dropped for simplicity, since it is a one di- 
mensional force. The details of this section can be found in 
Appendix D 1 of reference [ 1 ] . 



2. Equations of Motion for the Driven Parallel Plates 

The Unruh-Davies effect states that heating up of an ac- 
celerating conductor plate is proportional to its acceleration 
through the relation 



T 



2-Kck' ' 



where R is the plate acceleration. A one dimensional sys- 
tem of two parallel plates, shown in Figure ^2 can be used 
as a simple model to demonstrate the complicated sonolumi- 
nescense phenomenon for a bubble subject to a strong acoustic 
field. 

The dynamical force for the ID, linear coupled system can 
be expressed with equation ( 1641 . 



R\ — 771 i?i — 772-R2 — Crp, 

Ri - V3R2 - V1R1 = Zip, 



(65) 
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where 



Vl = m rp [Srp,2g a ,2 ~ Srp,3ffa,3j , 
Crp — Srp,39a 1 3m<rp Rdpr,cm,cti 

dip = —8i p ,ig a ,tm^ , 

Rl Rrp,crn,ai 
R2 Rip, cm, a? 



■ _ dR a " _ d 2 R a 

dt dt z 



(66) 



(67) 



Here i? 4 represents the center of mass position for the "Right 
Plate" and R2 represents the center of mass position for the 
"Left Plate" as illustrated in Figure m rp and mi p are the 
corresponding masses for the right and left plates. s rp .i and 
sjp i are the appropriate (±)signs that needs to be determined. 
For the mean time, the exact signs for s rp ^ and si P: i are not 
relevant, hence it is left as is. With a slight modification, equa- 
tion d65i for this linear coupled system can be written in the 
matrix form, (see Appendix D2 of reference 1 1]): 



and 















Ri = 


■/ 


R 3 dt', 


R2 = 


/ R 4 d?, 




J t 






'to 




' R 3 ~ 




m V2 




' R3' 


+ 


C,rp 






m m _ 




R4 




Zip _ 


/?„ 






Rt, 




1* 



(68) 



where 

Rl — i?3, i?2 = i?4, 

R3 =R\ = irp + mRi + V2R2 = £, rp + V1R3 + V2R4,, 
Ri = R 2 = 6 P + V3R2 + mRi = Zip + V3R4 + mR3- 

The matrix equation has the solutions (see Appendix D2 of 



reference fl]): 
R 



■rp,cm,a (^) ~~ 



A 4 (;t ) - m 0*o) 



dt' 



A3 (;*o) - vi (; £ o) 

ipn (t, t ) Rr P 

,cm,a (t ) + ^12 (t, t ) Rip .cm. a (to) 

exp([A 3 (;t ) + A 4 (;t )]t ) 
f* ^22 (?, t ) tj rp (?) ~ (?, t ) j lp (?) 

J to V>11 (f, to) V>22 (f , to) " V>12 (*', to) V>21 to) 
X Vll (t, t ) +^12 (t,t ) 

''' V-n (*', t ) 6 P (tQ - ^21 (f, t ) g rp (g) , 

t0 Vll (f, t ) ^22 (?, t Q ) ~ Vl2 (*', t ) V>21 (f, t ) ' 

(69) 



Rip, cm, a \v) 



A 4 (;to) - m (;to) 



.A3 (;t ) - rj! (;t ) 

l/j 2 l (t, t ) Rrp ,cm,a + V>22 (t, t ) % (to) 

exp([A 3 (;to) + A 4 (;to)]to) 
f* ^22 (*', t ) g rp (?) - jjy (?, t ) lip (?) 

J t V'll (?, to) ^22 (?, t ) ~ i>12 (?, t ) ihl (f, to) 
X 1p21 (t, t ) +"022 (Mo) 

''' V'll (*', t ) 6 P (tQ - ^21 (?, t ) t] rp (?) 

t0 ^11 (?, t ) 1p22 (?, t Q ) ~ Vl2 (f, t ) ^21 (f, t ) ' 

(70) 

where 



A3 = ^ h V 4 W ~ ? ' 3 J + 7727741 (71) 



?7i + »?3 



V'll (t, to) = i? 3 i?4 



1 2 

t Wi-»73] +?72?74, (72) 



A 4 - Tfl 



A3 - 771 
exp (A 4 t + A 3 t )] 



exp (A 3 t + A 4 t ) 



(73) 



V>12 (t, to) = R3R4 



V2 



A 3 - 771 



A 3 - 771 
exp (A 3 t + A 4 t 



exp (A 4 t + A 3 t ) 



(74) 



V"21 (t, to) = R3R4 

A 4 - 771 
772 

V>22 (t, to) = R3R 



A 4 - ??i 



exp (A 3 t + A 4 to) 



772 

exp (A 4 t + A 3 t ) 



(75) 



A 4 - 771 



A 3 - TJi 
exp (A 3 t + A 4 to)] 



exp (A 4 t + A 3 t ) 



(76) 
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The quantities R rp , C m,a and Ri PtCm ,a are the speed of the cen- 
ter of mass of "Right Plate" and the speed of the center of mass 
of the "Left Plate," respectively, and a defines the particular 
basis direction. The corresponding positions R rp , C m,a (t) and 
ftp, cm, a (t) are found by integrating equations d69i and J70i 
with respect to time, 



Edge of Universe 



Rr 



,a (*) 



A 4 (;to) - Vi (;t ) 



A3 (;*o) ~ Vi (;*o) 
ipn (r, t )R rp 

.cm .a (to) + 1P12 (r, t Q ) Ri p 

,cm,a (to) 

exp([A 3 (;to)+A 4 (;io)]^o) 

V>22 (f,to)Srp (t')-^X2 V,to)&p (?) 



I t0 Vll (*', t ) 1p22 (?, t ) ~ yJ!2 (f, t ) ^21 (*', t ) 

x tp n (t, t ) + tpi2 (t, t ) 

V'll (f, to) 6p (0 - V^21 (?, to) Zrp (?) 



dt' 



dt' 



/to ^11 (#, to) ^22 (f, t ) ~ ^12 (f, t ) ^21 (f, t ) 
X dr + Rrp,cm,a (to) , (77) 



ft 



lp,cm 



a(t) = 



A4 (;t ) - »?i (;t ) 



1 



A3 (;to) - m (;to) 

^21 (t, t )R rp .cm. a (to) + ^22 (r, t ) Rip .cm. a (to) 
exp([A 3 (;t )+A 4 (;to)]to) 

^22 (f.toKrp (O-fe (t',t )6 P (f) 



/to V'll (f , to) ^22 (f, t ) " V>12 (f , t ) ^21 (f, t ) 
X V21 (T) to) + V>22 (T, to) 

(tVtoKtg (tQ-jfex (t\tp)§rp(f) 



dt' 



dt' 



l t0 1/) U (f, to) lj>22 (?, to) ~ ^12 (f, t ) ^21 (*', t ) 

x dr + R lP:Cm , a (t ) • (78) 

The remaining integrations are straightforward and the ex- 
plicit forms will not be shown here. 

One may argue that for the static case, R rp ,cm,a (to) and 
Ri p ,cm,a (to) must be zero because the conductors seem to be 
fixed in position. This argument is flawed, for any wall to- 
tally fixed in position upon impact would require an infinite 
amount of energy. One has to consider the conservation of 
momentum simultaneously. The wall has to have moved by 
the amount AR wa u = R wa uAt, where At is the total dura- 
tion of impact, and R wa ii is calculated from the momentum 
conservation and it is non-zero. The same argument can be 
applied to the apparatus shown in Figure^ F° r that system 



llPuirtttai— P hoton || — —T~t ns sjf (to) 



and 



Rrp,cm,a (to) — ftp, 3 (to) + ftrp,2 (to) 
Rip .cm. a. (to) — ftrp,l (to) + ftp, 2 (to) 

For simplicity, assuming that the impact is always only in the 




Infinity 



Sphere vacuum— field Poynlinu vector field lines 



Figure 12: Boyer's configuration is such that a sphere is the only 
matter in the entire universe. His universe extends to the infinity, 
hence there are no boundaries. The sense of vacuum-field energy 
flow is along the radial vector f, which is defined with respect to the 
sphere center. 



normal direction, 



R 



rp,cm 



a {to 



and 



Rip .cm. a (to) 



mipC 



\ft'n e ,3 (to) - H'n B ,2 ( f o)\ 



\l~?n s ,l (to) - 7~L'n s ,2 (*o)|| > 



where the differences under the magnitude symbol imply field 
energies from different regions counteract the other. The de- 
tails of this section can be found in Appendix D2 of reference 

a. 



III. RESULTS AND OUTLOOK 

The results for the sign of Casimir force on non-planar ge- 
ometric configurations considered in this investigation will 
eventually be compared with the classic repulsive result ob- 
tained by Boyer decades earlier. For this reason, it is worth 
reviewing Boyer's original configuration as shown in Figure 

T. H. Boyer in 1968 obtained a repulsive Casimir force re- 
sult for his charge-neutral, hollow spherical shell of a perfect 
conductor |12|. For simplicity, his sphere is the only object 
in the entire universe and, therefore, no external boundaries 
such as laboratory walls, etc., were defined in his problem. 
Furthermore, the zero-point energy flow is always perpendic- 
ular to his sphere. Such restriction constitutes a very stringent 
condition for the material property that a sphere has to meet. 
For example, if one were to look at Boyer's sphere, he would 
not see the whole sphere; but instead, he would see a small 
spot on the surface of a sphere that happens to be in his line of 
sight. This happens because the sphere in Boyer's configura- 
tion can only radiate in a direction normal to the surface. One 
could equivalently argue that Boyer's sphere only responds to 
the approaching radiation at normal angles of incidence with 
respect to the surface of the sphere. When the Casimir force 
is computed for such restricted radiation energy flow, the re- 
sult is repulsive. In Boyer's picture, this may be attributed to 
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Figure 13: Manufactured sphere, in which two hemispheres 
are brought together, results in small non-spherically symmetric 
vacuum-field radiation inside the cavity due to the configuration 
change. For the hemispheres made of Boyer's material, these fields 
in the resonator will eventually get absorbed by the conductor result- 
ing in heating of the hemispheres. 



the fact that closer to the origin of a sphere, the spherically 
symmetric radiation energy flow becomes more dense due to 
the inverse length dependence, and this density decreases as 
it gets further away from the sphere center. This argument, 
however, seems to be flawed because it inherently implies ex- 
istence of the preferred origin for the vacuum fields. As an 
illustration, Boyer's sphere is shown in Figure^l F° r the rest 
of this investigation, "Boyer's sphere" would be strictly re- 
ferred to as the sphere made of material with such a property 
that it only radiates or responds to vacuum-field radiations at 
normal angle of incidence with respect to its surface. 

The formation of a sphere by bringing together two nearby 
hemispheres satisfying the material property of Boyer's 
sphere is illustrated in Figure[0] Since Boyer's material prop- 
erty only allow radiation in the normal direction to its surface, 
the radiation associated with each hemisphere would neces- 
sarily go through the corresponding hemisphere centers. For 
clarity, let us define the unit radial basis vector associated with 
the left and right hemispheres by t\l and ?r, respectively. If 
the hemispheres are made of normal conductors the radiation 
from one hemisphere entering the other hemisphere cavity 
would go through a complex series of reflections before es- 
caping the cavity. Here, a conductor with Boyer's stringent 
material property is not considered normal. Conductors that 
are normal also radiate in directions non-normal to their sur- 
face, whereas Boyer's conductor can only radiate normal to 
its surface. Due to the fact that Boyer's conducting materials 
can only respond to radiation impinging at a normal angle of 
incidence with respect to its surface, all of the incoming radi- 
ation at oblique angles of incidence with respect to the local 
surface normal is absorbed by the host hemisphere. This sug- 
gests that for the hemisphere-hemisphere arrangement made 
of Boyer's material shown in Figure ^] the temperature of 
the two hemispheres would rise indefinitely over time. This 
does not happen with ordinary conductors. This suggests that 
Boyer's conducting material, of which his sphere is made, is 
completely hypothetical. Precisely because of this material 
assumption, Boyer's Casimir force is repulsive. 

For the moment, let us relax the stringent Boyer's mate- 
rial property for the hemispheres to that of ordinary conduc- 
tors. For the hemispheres made of ordinary conducting ma- 



Figure 14: The process in which a configuration change from 
hemisphere-hemisphere to sphere inducing virtual photon in the di- 
rection other than r is shown. The virtual photon here is referred to 
as the quanta of energy associated with the zero-point radiation. 



terials, there would result a series of reflections in one hemi- 
sphere cavity due to those radiations entering the cavity from 
nearby hemisphere. For simplicity, the ordinary conducting 
material referred to here is that of perfect conductors with- 
out Boyer's hypothetical material property requirement. Fur- 
thermore, only the radiation emanating normally with respect 
to its surface is considered. The idea is to illustrate that the 
"normally emanated radiation" from one hemisphere results 
in elaboration of the effects of "obliquely emanated radiation" 
on another hemisphere cavity. Here the obliquely emanated 
radiation means those radiation emanating from a surface not 
along the local normal of the surface. 

When two such hemispheres are brought together to form a 
sphere, there would exist some radiation trapped in the sphere 
of which the radiation energy flow lines are not spherically 
symmetric with respect to the sphere center. To see how a 
mere change in configuration invokes such non-spherically 
symmetric energy flow, consider the illustration shown in Fig- 
ure^J For clarity, only one "normally emanated radiation" 
energy flow line from the left hemisphere is shown. When one 
brings together the two hemispheres just in time before that 
quantum of energy escapes the hemisphere cavity to the right, 
the trapped energy quantum would continuously go through 
series of complex reflections in the cavity obeying the reflec- 
tion law. But how fast or how slow one brings in two hemi- 
spheres is irrelevant in invoking such non-spherically sym- 
metric energy flow because the gap 8 can be chosen arbitrarily. 
Therefore, there would always be a stream of energy quanta 
crossing the hemisphere opening with £ ^ as shown in Fig- 
ure [3] In other words, there is always a time interval At 
within which the hemispheres are separated by an amount 5 
before closure. The quanta of vacuum-field radiation energy 
created within that time interval At would always be satis- 
fying the condition £ ^ 0, and this results in reflections at 
oblique angle of incidence with respect to the local normal of 
the walls of inner sphere cavity. Only when the two hemi- 
spheres are finally closed, would then £ = and the radiation 
energy produced in the sphere after that moment would be 
spherically symmetric and the reflections would be normal to 
the surface. However, those trapped quantum of energy that 
were produced prior to the closure of the two hemispheres 
would always be reflecting from the inner sphere surface at 
oblique angles of incidence. 
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Poynting vector field lines originating 
from lab boundaries 



■ Lab boundary _ 



Figure 15: A realistic laboratory has boundaries, e.g., walls. These 
boundaries have effect similar to the field modes between two paral- 
lel plates. In 3D, the effects are similar to that of a cubical laboratory, 
etc. 
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Figure 16: The schematic of sphere manufacturing process in a real- 
istic laboratory. 



Unlike Boyer's ideal laboratory, realistic laboratories have 
boundaries made of ordinary material as illustrated in Figure 
IT31 One must then take into account, when calculating the 
Casimir force, the vacuum-field radiation pressure contribu- 
tions from the involved conductors, as well as those contribu- 
tions from the boundaries such as laboratory walls, etc. We 
will examine the physics of placing two hemispheres inside 
the laboratory. 

For simplicity, the boundaries of the laboratory as shown 
in Figure^]are assumed to be simple cubical. Normally, the 
dimension of conductors considered in Casimir force experi- 
ment is in the ranges of microns. When this is compared with 
the size of the laboratory boundaries such as the walls, the 
walls of the laboratory can be treated as a set of infinite paral- 
lel plates and the vacuum-fields inside the the laboratory can 
be treated as simple plane waves with impunity. 

The presence of laboratory boundaries induce reflection of 
energy flow similar to that between the two parallel plate ar- 
rangement. When the two hemisphere arrangement shown in 
Figure ^] is placed in such a laboratory, the result is to elab- 
orate the radiation pressure contributions from obliquely inci- 
dent radiations on external surfaces of the two hemispheres. If 
the two hemispheres are made of conducting material satisfy- 
ing Boyer's material property, the vacuum-field radiation im- 
pinging on hemisphere surfaces at oblique angles of incidence 
would cause heating of the hemispheres. It means that Boyer's 
hemispheres placed in a realistic laboratory would continue to 
rise in temperature as a function of time. However, this does 
not happen with ordinary conductors. 

If the two hemispheres are made of ordinary perfect con- 
ducting materials, the reflections of radiation at oblique angles 
of incidence from the laboratory boundaries would elaborate 



on the radiation pressure acting on the external surfaces of two 
hemispheres at oblique angles of incidence. Because Boyer's 
sphere only radiates in the normal direction to its surface, or 
only responds to impinging radiation at normal incidence with 
respect to the sphere surface, the extra vacuum-field radiation 
pressures considered here, i.e., the ones involving oblique an- 
gles of incidence, are missing in his Casimir force calculation 
for the sphere. 



A. Results 

T. H. Boyer in 1968 have shown that for a charge-neutral, 
perfect conductor of hollow spherical shell, the sign of the 
Casimir force is positive, which means the force is repulsive. 
He reached this conclusion by assuming that all vacuum-field 
radiation energy flows for his sphere are spherically symmet- 
ric with respect to its center. In other words, only the wave 
vectors that are perpendicular to his sphere surface were in- 
cluded in the Casimir force calculation. In the following sec- 
tions, the non-perpendicular wave vector contributions to the 
Casimir force that were not accounted for in Boyer's work are 
considered. 



1. Hollow Spherical Shell 

As shown in Figure [^] the vacuum-field radiation imparts 
upon a differential patch of an area dA on the inner wall of the 
conducting spherical cavity a net momentum of the amount 



^Pinner ^H/\k inner R s,l) R s,0^ 



2nnh cos 9 r 



R s , 2 lrLA> 



' Ji s,2 



R s ,i 



R 1 



8,1, 



where 



o < e inc < tt/2, 

n= 1,2,3, ••• . 

Here Ak 1 inner (l R's.i, R's.tjj is from equation 04l >. The an- 
gle of incidence 0j nc is from equation i2U : R sA (r^, A' s 

and i? Sj 2 (r'iiA's follow the generic form shown in equa- 
tion d35l . 

Similarly, the vacuum-field radiation imparts upon a differ- 
ential patch of an area dA on the outer surface of the conduct- 
ing spherical shell a net momentum of the amount 

Ap ou ter — — HAk r outer 



2 

= -2ft 



(:,R' sA + aR' sA 

COS Oi nc R' s i, 



where 



o < e mc < ?r/2, 
n= 1,2,3,- •• . 
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Figure 17: The vacuum-field wave vectors k'i.b and k'ij impart a 
net momentum of the magnitude ||pnet|| = h |fc*'i,t> — /2 on 

differential patch of an area dA on a conducting spherical surface. 



aR' Si i ) is from equation (|36 



The net average force per unit time, per initial wave vector 
direction, acting on differential element patch of an area dA is 
given by 



•F, s,avq lim 



Ap 



At 



At 



or 



s,avg 



hf 



4,2 (r^A^)-^^,!^) 
x 2h cos 9 inc R' s i, 



where 

o < e inc < tt/2, 

u= 1,2,3,- •• . 

Notice that T Syavg is called a force per initial wave vector di- 
rection because it is computed for k'i,b an d k'ij along specific 
initial directions. Here k'ij denotes a particular initial wave 
vector k'i entering the resonator at R St o as shown in Figure[5] 
The subscript b for k'ij denotes the bounded space inside the 
resonator. The k'ij denotes a particular initial wave vector k'i 
impinging upon the surface of the unbounded region of sphere 
at point R' s .i + aR' S: i as shown in Figure|5] The subscript / 
for k'ij denotes the free space external to the resonator. 
Because the wave vector k'ij resides in free (unbounded) 



space, its magnitude 



hf 



can take on a continuum of 



allowed modes. The wave vector fc', h however resides in 



bounded region, hence 



is restricted by the relation 



L 



R 



,, 2 (^,a; 2 ) -iLtfaKi) 



The free space limit is the case where the radius of the sphere 



becomes very large. Therefore, by designating k'ij as 



ij 



lim 



rnr 



R 



s,2 



(r;,A\ 2 )-4,i(r<,A' s ,i) 



and summing over all allowed modes, the total average force 
per unit time, per initial wave vector direction, per unit area is 
given by 



*F : s,avq — R s.l 



E 

n=l 



t2H COS 6ir, 



— Jim 



4,2 (r-U' s , 2 )-4,i ( r 'nKi) 
m:2hcos #,„ r 



" \\R.,2 (r'i, A' 



s.2 



Rs.i 



In the limit r[ — > oo, the second summation to the right can 
be replaced by an integration, YlnLi ~* Io° dn. Hence, we 
have 



•F s,avq — R s.l 



E 



2Hmr cos 6i„ c 



n=l \\Rs.2 



lim 



faA'^-R^ (riA' s>1 ) 
2hmr cos 9i nc dn 



4,2 (r<,A*' s , 2 )-4,i (<,K,i) 



or with the following substitutions, 



k' 



hf 



RsA K.A'i 



dn 



R s 



{r'iX^-Rs,!^,^ 



the total average force per unit time, per initial wave vector 
direction, per unit area is written as 



■Ts.avg — 7 II ^ / \ -> / -.. \ 

[^i\\Rs,2 (^A^)-i? Sjl (r^A' Sil ) 
- ijin^ (1^,2 (^,A' s 2 ) - R Sil 

Kjdklj^j 2hcos9 inc R' sA , (79) 



where < 6i nc < n/2 and n = 1, 2, 3, ■ ■ • . The total aver- 
age vacuum-field radiation force per unit time acting on the 
uncharged conducting spherical shell is therefore 



F, 



s.total 



E^ 



s,avg ' dS sphere 
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F, 



s. total — 



_E^ 

{k'i,b,k' i,/,R' 5 ,o} 



2nnh cos 6;. 



k i jdk i j — cos( 



^ s,l ' dS sphere^ 

(80) 



where dS sp here is a differential surface element of a sphere 
and the integration J s is over the spherical surface. The term 
R's.o is the initial crossing point inside the sphere as defined in 
equation @. The notation Ylw i & f R' } ' m P^y me sum- 
mation over all initial wave vector directions for both inside 
(js'i.bj and outside (k'ij^j of the sphere, over all crossing 

points given by i?' s o- 

It is easy to see that !F s ,avg of equation J79l > is an "unregu- 
larized" ID Casimir force expression for the parallel plates 
(see the vacuum pressure approach by Milonni, Cook and 
Goggin 1 22]). It becomes more apparent with the substitution 
At = d/c. An application of the Euler-Maclaurin summation 
formula 1 26, 27] leads to the regularized, finite force expres- 
sion. The force T s ,avg is attractive because 

cos 8 mn > 



and 



E 



n=l \\Rs, 2 (r' v A' 
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Rs,i (r>,K< s l ) 
- Jim ( R^UJ'^-R^U,^ 



where 



R 



8,2 



.A 



8,2 



- R 



A 



s.l 



is a constant for a 



given initial wave k'i t b and the initial crossing point R' s .o in 
the cross-section of a sphere (or hemisphere). The total av- 
erage force Fs^otai, which is really the sum of T s . aV g over 
all R's.o and all initial wave directions, is therefore also at- 
tractive. For the sphere configuration of Figure |5J where the 
energy flow direction is not restricted to the direction of local 
surface normal, the Casimir force problem becomes an exten- 
sion of infinite set of parallel plates of a unit area. 



2. Hemisphere-Hemisphere and Plate-Hemisphere 

For the hemisphere-hemisphere and plate-hemisphere con- 
figurations, the expression for the total average force per unit 



time, per initial wave vector direction, per unit area is identical 
to that of the hollow spherical shell with modifications, 



•Fh.avq — R h. 
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i \\R 
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x / kifdki f 
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2fi,cos 6j 



(81) 



where 9i nc < it/2 and n = 1,2,3, •■ • . The inci- 
dence angle 6. mc is from equation fl21i : Rh,i {r'i, A^ and 

Rh,2 (r'i, AJj 2 ^ follow the generic form shown in equation 
j39t . This force is attractive for the same reasons as dis- 
cussed previously for the hollow spherical shell case. The 
total radiation force averaged over unit time, over all possible 
initial wave vector directions, acting on the uncharged con- 
ducting hemisphere-hemisphere (plate-hemisphere) surface is 
given by 



F, 



h, total 



{k' i,b-k' ij ,R' h l0 } 
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2nnh cos Qj. 



i K.At a )-R h , x r^,A^ 
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' it ,rt 2h 

< I k i jdk i j — cos( 
lo ' 7r 



R'h,l ' dSfiemisphere, 

(82) 



where dShemisphere is now a differential surface element of 
a hemisphere and the integration J s is over the surface of the 

hemisphere. The term R'h.o is the initial crossing point of the 
hemisphere opening as defined in equation (0}. The notation 
Ylw b i? f R' h o} i m Ply me summation over all initial wave 

vector directions for both inside (k'i.bj and outside (k'ijj 
of the hemisphere-hemisphere (or the plate-hemisphere) res- 
onator, over all crossing points given by R'h,o- 

It should be remarked that for the plate-hemisphere config- 
uration, the total average radiation force remains identical to 
that of the hemisphere-hemisphere configuration only for the 
case where the gap distance between plate and the center of 
hemisphere is more than the hemisphere radius r[. When the 

plate is placed closer, the boundary quantization length 
must be chosen carefully to be either 
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They are illustrated in Figure [9] The proper one to use is the 

smaller of the two. Here R p (r^, Ay J is from equation ( 149 \ 

of Appendix C3 and Nh.max is defined in equation (I40i of 
Appendix C2. 

B. Interpretation of the Result 

Because only the specification of boundary is needed in 
Casimir's vacuum-field approach as opposed to the use of 
a polarizability parameter in Casimir-Polder interaction pic- 
ture, the Casimir force is sometimes regarded as a configu- 
rational force. On the other hand, the Casimir effect can be 
thought of as a macroscopic manifestation of the retarded van 
der Waals interaction. And the Casimir force can be equiva- 
lently approximated by a summation of the constituent molec- 
ular forces employing Casimir-Polder interaction. This prac- 
tice inherently relies on the material properties of the involved 
conductors through the use of polarizability parameters. In 
this respect, the Casimir force can be regarded as a material 
dependent force. 

Boyer's material property is such that the atoms in his con- 
ducting sphere are arranged in such manner to respond only to 
the impinging radiation at local normal angle of incidence to 
the sphere surface, and they also radiate only along the direc- 
tion of local normal to its surface. When the Casimir force is 
calculated for a sphere made of Boyer's fictitious material, the 
force is repulsive. Also, in Boyer's original work, the labora- 
tory boundary did not exist. When Boyer's sphere is placed 
in a realistic laboratory, the net Casimir force acting on his 
sphere becomes attractive because the majority of the radia- 
tion from the laboratory boundaries acts to apply inward pres- 
sure on the external surface of sphere when the angle of inci- 
dence is oblique with respect to the local normal. If the sphere 
is made of ordinary perfect conductors, the impinging radia- 
tion at oblique angles of incidence would be reflected. In such 
cases the total radiation pressure applied to the external local- 
sphere-surface is twice the pressure exerted by the incident 
wave, which is the force found in equation ( 18 Ot of the previ- 
ous section. However, Boyer's sphere cannot radiate along the 
direction that is not normal to the local-sphere-surface. There- 
fore, the total pressure applied to Boyer's sphere is half of the 
force given in equation ( 1801 1 of the previous section. This pe- 
culiar incapability of emission of a Boyer's sphere would lead 
to the absorption of the energy and would cause a rise in the 
temperature for the sphere. Nonetheless, the extra pressure 
due to the waves of oblique angle of incidence is large enough 
to change the Casimir force for Boyer's sphere from being 
repulsive to attractive. The presence of the laboratory bound- 
aries only act to enhance the attractive aspect of the Casimir 
force on a sphere. The fact that Boyer's sphere cannot irradi- 
ate along the direction that is not normal to the local-sphere- 
surface, whereas ordinary perfect conductors irradiate in all 
directions, implies that his sphere is made of extraordinarily 
hypothetical material, and this may be the reason why the re- 
pulsive Casimir force have not been experimentally observed 
to date. 




Figure 18: To deflect away as much possible the vacuum-field ra- 
diation emanating from the laboratory boundaries, the walls, floor 
and ceiling are constructed with some optimal curvature to be deter- 
mined. The apparatus is then placed within the "Apparatus Region." 



In conclusion, (1) the Casimir force is both boundary and 
material property dependent. The particular shape of the con- 
ductor, e.g. sphere, only introduces the preferred direction 
for radiation. For example, radiations in direction normal to 
the local surface has bigger magnitude than those radiating 
in other directions. This preference for the direction of ra- 
diation is intrinsically connected to the preferred directions 
for the lattice vibrations. And, the characteristic of lattice 
vibrations is intrinsically connected to the property of ma- 
terial. (2) Boyer's sphere is made of extraordinary conduct- 
ing material, which is why his Casimir force is repulsive. (3) 
When the radiation pressures of all angles of incidence are in- 
cluded in the Casimir force calculation, the force is attractive 
for charge-neutral sphere made of ordinary perfect conduc- 
tor. And, lastly, (4) the Casimir force problem involving any 
non-planar geometric boundary configurations can always be 
reduced to a series of which the individual terms in the series 
are that of the parallel plates problem. Because each terms 
in the series are that of parallel plates problem, the summa- 
tion over these individual terms result in an attractive Casimir 
force. 



C. Suggestions on the Detection of Repulsive Casimir Force 
for a Sphere 

The first step in detecting the repulsive Casimir force for 
a spherical configuration is to find a conducting material 
that most closely resembles the Boyer's material to construct 
two hemispheres. It has been discussed previously that even 
Boyer's sphere can produce attractive Casimir force when the 
radiation pressures due to oblique incidence waves are in- 
cluded in the calculation. Therefore, the geometry of the lab- 
oratory boundaries have to be chosen to deflect away as much 
as possible the oblique incident wave as illustrated in Figure 
ITSl Once these conditions are met, the experiment can be con- 
ducted in the region labeled "Apparatus Region" to observe 
Boyer's repulsive force. 
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Figure 19: The original bubble shape shown in dotted lines and the 
deformed bubble in solid line under strong acoustic field. 



D. Outlook 



Figure 20: The vacuum-field radiation energy flows are shown for 
closed and unclosed hemispheres. For the hemispheres made of 
Boyer's material, the non-radial wave would be absorbed by the 
hemispheres. 



The Casimir effect has influence in broad range of physics. 
Here, we list one such phenomenon known as "sonolumi- 
nescense," and, finally conclude with the Casimir oscillator. 



1. Sonoluminescense 

The phenomenon of sonoluminescense remains a poorly 
understood subject to date |28, 29]. When a small air bubble 
of radius ~ 10~ 3 cm is injected into water and subjected to 
a strong acoustic field of ~ 20 kHz under pressure roughly 
~ 1 atm, the bubble emits an intense flash of light in the 
optical range, with total energy of roughly ~ 10 7 eV. This 
emission of light occurs at minimum bubble radius of roughly 
~ 10 ~ 4 cm. The flash duration has been determined to be on 
the order of 100 ps 1301 l3ll 13211 . It is to be emphasized that 
small amounts of noble gases are necessary in the bubble for 
sonoluminescense. 

The bubble in sonoluminescense experiment can be 
thought of as a deformed sphere under strong acoustic 
pressure. The dynamical Casimir effect arises due to the 
deformation of the shape; therefore, introducing a modifica- 



tion to L21 = 



i?2 — Ri 



from that of the original bubble 



shape. Here L21 is the path length for the reflecting wave 
in the original bubble shape. In general L21 = £21 (t) — 
\R2 (n (t) ,e(t),<t> (t)) - Ri (n (t) ,6(t),<f> (i))|| . From 
the relations found in this work for the reflection points 
Ri (n (t) , (t) , <j} (*)) and R N (r ( (t) , 6 (t) , <f> (*)) , together 
with the dynamical Casimir force expression of equation 
(1571 . the amount of initial radiation energy converted into 
heat energy during the deformation process can be found. 
The bubble deformation process shown in FigurefH^is a three 
dimensional heat generation problem. Current investigation 
seeks to determine if the temperature can be raised sufficiently 
to cause deuterium-tritium (d-t) fusion to occur, which could 
provide an alternative approach to achieve energy generation 
by this d-t reaction (threshold ~ 17 KeV) ll33ll . Its theoretical 
treatment is similar to that discussed on the ID problem 
shown in Figure ITOl 



2. Casimir Oscillator 

If one can create a laboratory as shown in Figure and 
place in the laboratory hemispheres made of Boyer's mate- 
rial, then the hemisphere-hemisphere system will execute an 
oscillatory motion. When two such hemispheres are sepa- 
rated, the allowed wave modes in the hemisphere-hemisphere 
confinement would no longer follow Boyer's spherical Bessel 
function restriction. Instead it will be strictly constrained 



by the functional relation of 



i?2 — Ri 



where Ri and R2 

are two neighboring reflection points. Only when the two 
hemispheres are closed, would the allowed wave modes obey 
Boyer's spherical Bessel function restriction. 

Assuming that hemispheres are made of Boyer's material 
and the laboratory environment is that shown in Figure^] the 
two closed hemispheres would be repulsing because Boyer's 
Casimir force is repulsive. Once the two hemispheres are sep- 
arated, the allowed wave modes are governed by the inter- 
nal reflections at oblique angle of incidence. Since the hemi- 
spheres made of Boyer's material are "infinitely unrespon- 
sive" to oblique incidence waves, all these temporary non- 
spherical symmetric waves would be absorbed by the Boyer's 
hemispheres and the hemispheres would heat up. The two 
hemispheres would then attract each other and the oscillation 
cycle repeats. Such a mechanical system may have applica- 
tion. 
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